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Abstract 

All mutually unbiased bases in dimension six consisting of product states 
only are constructed. Several continuous families of pairs and two triples of 
mutually unbiased product bases are found to exist but no quadruple. The 
exhaustive classification leads to a proof that a complete set of seven mutu- 
ally unbiased bases, if it exists, cannot contain a triple of mutually unbiased 
product bases. 



1 Introduction 

Mutually unbiased (MU) bases [iHS] have attracted interest in recent years because 
their properties seem to depend dramatically on the dimension d of the quantum 
system in hand. If the dimension is given by a prime number the state space 
is known to accommodate a complete set of (p + 1) MU bases. Each of these 
bases consists of p orthonormal states \ ja) G C^, with constant overlap of 1/p across 
different bases, 

\{ia\h)? = -{'^ - 5ab) + Sjk^ab, j , k = . . . p - 1 , a, b = . . . p . (1) 

p 

Complete sets of MU bases also exist for quantum systems with dimension d = 
p", where n is a positive integer. However, for "composite" dimensions such as 
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d = did2 G {6, 10, 12, . . .} complete sets of MU bases seem to be absent. In spite of 
considerable numerical searches P|[5] , computer-algebraic efforts |6|I7] , and numerical 
calculations with rigorous error bounds, only three MU bases have been found in 
dimension six, four less than the maximally allowed number |8]. Thus, the six- 
dimensional state space of a qubit-qutrit system appears to differ structurally from 
the state space of a pair of qubits {d = 4) or a pair of qutrits [d = 9). 

One of the few known results in dimension d = 6 is the impossibility to extend, by 
more that one further MU basis, the pair of MU bases consisting of the standard 
basis and its dual, the Fourier basis [5]. Thus, triples of MU bases are the largest 
sets to be found in this way. Another, more recent result [8] states that the Fourier 
family of Hadamard matrices together with the identity cannot be extended to a MU 
quadruple. These initial pairs, after non-local equivalence transformations, consist 
of product states only, a fact which has received little attention. 

Upon reflection, it seems worthwhile to systematically study MU bases in composite 
dimensions which contain only product states. In the present paper we carry out a 
comprehensive study of MU product bases in dimensions six, complementing studies 
devoted to the entanglement structure of complete sets of MU bases [9]-[TT]. 

More specifically, we will derive an exhaustive list of MU product bases in dimension 
six. The restriction to product states goes hand in hand with local equivalence 
transformations, or LETs, consisting of local (anti-) unitary transformations. We 
will find that in the space C^, there is a considerable number of inequivalent 

product bases, a limited set of families of MU product pairs and just two triples of 
MU product bases. No larger MU product constellations exist. This result effectively 
limits the number of MU product bases contained in a hypothetical complete set of 
MU bases in dimension six. 

The argument will unfold as follows. In Sec. 2 we introduce MU product bases, 
specify all local (anti-) unitary transformations which map a given set of MU product 
states to an equivalent one, and summarise relevant properties of MU bases in 
dimensions two and three. Then, in Sec. 3, we derive all inequivalent product 
bases in and C^. Sec. 4 has two results on product vectors required to be MU to 
certain given sets of MU product vectors. These results will be important tools to 
enumerate all pairs and triples of MU bases in dimension four (Sec. 5) and dimension 
six (Sec. 6). This classification allows us to conclude, as shown in Sec. 7, that no 
MU product triple can be part of a complete set of seven MU bases in d = Q. The 
final section summarises our findings. 

Readers mainly interested in the results relevant to dimension six are advised to 
immediately proceed to Sec. 6 after having familiarized themselves with the concept 
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of mutually unbiased product bases presented in Sec. 2. 



2 MU product bases 

From now on, we will consider quantum systems consisting of two subsystems with 
prime dimensions p and q, where p < q. The state space of such a bipartite system is 
given by the Hilbert space ® of dimension d = pq. Since g is a prime number, 
there is a complete set of MU bases of C^, and we will denote its q{q + 1) states by 

|Jb)eO, J = 0...g-1, 6 = 0...g. (2) 

The q states {| J;,)} form one orthonormal basis labelled by b, and states taken from 
two distinct bases are MU, in analogy to Eq. ([T]). Given complete sets of MU bases 
in and C^, respectively, we now construct (p+ 1) MU product bases of the space 
(g) C^. To do so, we pair each MU basis of the space C with a (different) basis 
of and, within each pair, we tensor each state of the first basis with a (different) 
state of the second one. This procedure results in pq{p + 1) product states 

\ja)(^\Ja) = \ja,Ja) , (3) 

forming (p+ 1) MU bases {\ja, Jo), a = . . .p} of the space C^C^C. This is evident 
upon calculating the overlaps 

IO.,J#.,i^.)r = |(j#,)ri(J.|i^,)r=| i^^^^ = (4) 

which are the conditions for bases to be MU in a space of dimension pq. 

One can construct MU product bases of the type given in Eq. ([3]) using Heisenberg- 
Weyl (HW) operators. In dimension p, with p prime, the HW cyclic shift (modulo 
p) and phase operators Xp and Zp, respectively, are defined as 

Xp\j) = \j + 1) and Zp\j)=u^\j), (5) 

where u = e^'^'^lf is a p^"^ root of unity and is the standard basis with j = 

. . .p — 1. Since p is prime, one can construct a complete set of (p + 1) MU bases 
from the eigenbases of the operators Xp{Z^^ for < £ < p — 1 |12| . 

For the composite dimension d = pq, we can build a set of {p+ 1) MU product bases 
of the Hilbert space ® C'' with the operators Xp and Zp acting on the space C^, 
and Xq and Zg on the space C^. For example, the eigenbases of the operators Xp^Xg 
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and ® Zq form two MU product bases, which we call a Heisenberg-Weyl pair. 
One can also construct HW bases with the operators X^^ and on the space C^', 
however, these do not necessarily form product bases. Since we are concerned with 
product bases in this paper, we define the HW operators on the space ® such 
that their eigenstates are product states. Note that we do not limit the construction 
of MU bases to the eigenbases of HW operators, i.e. {|ja;<^a)} in © can be any 
product basis. 

Each basis {|ja,</a)} is a direct product basis of the space ® C"^ since each 
state \3f^^3 = . . .p — 1, of the a^^ basis in is multiplied with every state 
\Ja), J = ... g — 1, of the a**^ basis of C^. Direct product bases are, however, only a 
subset among all product bases: indirect product bases jl^ result if the states being 
tensored stem from more than one basis of the space (or C^). The four states 

{|0„0,),|0„1,),|1„0,),|U,1,)} (6) 

provide a simple example of an indirect product basis in dimension four since two 
different bases of the second space, and occur in the construction. 

The matrix representation of a direct product basis in dimension d = pq is given by 
the tensor product of two matrices, each representing a basis of the spaces C and 
C"^, respectively. The matrix representation of an indirect product basis cannot be 
written as a tensor product of two matrices. 

Conceptually, the distinction between direct and indirect product bases is not linked 
to MU bases: instead of using and in (E]) any other pair of bases of 

would also define an indirect product basis. Indirect product bases are important 
since they have been found to exhibit a degree of non-locality in the absence of 
entanglement |13| . 

In this paper, we will be concerned exclusively with product bases of the spaces 
and (8) C^. To simplify the construction of all different MU product 
bases, we will now introduce equivalence relations which respect the structure of 
product states, followed by a brief reminder of the properties of MU bases in and 
following conventions used in |14] . 

2.1 Local equivalence transformations 

Given a set of MU bases on the space C^, we obtain another set by applying one sin- 
gle unitary transformation to all states simultaneously. The scalar products between 
the states of the MU bases do not change under this transformation so that we deal 
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indeed with a second set of MU bases, factually different from the initial set but 
equivalent to it. By not distinguishing between equivalent MU bases, their enumer- 
ation is greatly simplified. When representing MU bases by Hadamard matrices, the 
concept of a standard (or dephased) form emerges naturally (see [15], for example). 
To enumerate all MU product bases it will be helpful not to distinguish those sets 
of MU product bases which can be transformed into each other by local equivalence 
transformations, or LETs, for short. LETs are defined by the requirement that they 
preserve the product structure of all states. If there is no LET transforming two 
given sets of MU product states into each other they will be called locally inequiva- 
lent, or just inequivalent. It may still be possible to transform them into each other 
by non-local transformations. 

We now list all LETs for a bipartite quantum system with Hilbert space C ® C^. 
Suppose we are given sets of (r + 1) MU bases {Bq, Bi, . . . , Br} that contain only 
product states. Explicitly, the p"^ basis, with p = . . . r, consists oi d = pq product 
states \np,Np), n = N e {l,2,...,d}, where \np) G C*' and \Np) G O. Any 
combination of the following five operations maps the given set of MU bases into a 
locally equivalent set: 

1. a local unitary transformation u®U effecting 

Bp^B'^ = u<^UBp = [..., \unp) ® |f/iV,), . . .} , (7) 

which leaves invariant the value of all scalar products; 

2. the multiplication of all states within a basis by possibly different phase factors 
such that 

S,^S;=|...,e^*"|n„iV,),...}; (8) 

these transformations exploit the fact that the overall phase of a quantum state 
has no physical significance and automatically drops out from the conditions 
defining MU bases. It is worth noting that a single phase factor e*''^ can dephase 
both states of a product: let = 0' + 0" to find e^'^\np,Np) = {e'^'^'\np)) ® 
{e^^"\Np)); 

3. permutations of the product states within each basis; as an example, consider 
the permutation of states \np,Np) and \n'p,Np) in the p^^ basis 

[...,\np,Np),...,\n'p,N'^),...] {. . . , l^^, iV;), . . . , |n„ iV,), . . .} , (9) 
which amounts to relabelling the elements within each basis; 
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4. the local complex conjugations k®! and i®K (anti-unitary operations defined 
with respect to the standard bases in and C, respectively), and thus their 
product k ® K] for example, applying k ® I 

Sp^s; = {...,K,iv,),...}, (10) 

swaps all scalar products resulting from the first factors without changing their 
numerical values; 

5. pairwise exchanges of two bases, which amounts to relabelling the bases. 

We now briefly discuss some important properties of LETs. First, they represent a 
true subset of all equivalence transformations in a space of dimension pq: no LET 
maps an indirect product basis to a direct one while a general unitary equivalence 
transformation can send any orthonormal basis to any other. Second, we will find 
indirect product bases which cannot be transformed into each other by LETs, i.e. 
locally inequivalent product bases. As a result, the idea of a unique standard or 
dephased form of MU bases is less straightforward for MU product bases. We define 
a standard form in the following way: the first basis Bq, be it direct or indirect, 
contains the states {Ijz)} of the space and the states {| J^)} of the space O; the 
second basis Bi contains the state {O^, Ox), and all other states in the remaining bases 
are dephased using the transformation defined in ([H]). Superficially, LETs remind 
one of local operations with classical communication, or LOCCs jlS]. However, the 
presence of anti-unitary operations rather suggests a link with Wigner's theorem 
about symmetry transformations leaving transition probabilities invariant [T7J, for 
the special case of a universe populated with product states only. Finally, it is 
straightforward to generalise LETs to ri-partite systems residing in product states 
only. 

It is often convenient to represent an MU product basis in C^^ as a complex Hadamard 
matrix of dimension {pqxpq), with each product state corresponding to one column. 
The bases {Bq, Bi, . . . , Br} then turn into a set of (r -|- 1) matrices, on which the five 
transformations above act in the following way. The first LET is a local unitary, 
given by the Kronecker product of two unitary matrices, applied to all matrices 
from the left; the second LET corresponds to diagonal unitary transformations act- 
ing from the right; unitary permutation matrices acting from the right implement 
the third type of LET, while the effect of the local complex conjugations must be 
worked out by writing down each product state individually. 
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2.2 MU bases in dimensions two and three 



Given a pair of MU bases in the vector space C^, we can always map the first basis 
to the standard basis by a suitable unitary transformation u G SU(2). Being 

MU to the first basis, the states of the second basis now must have the form 

l«) = ^(|0.)+e^'|l.)) = r,|+), \a^)=f,\-), (11) 

where {|±)} = is the x-eigenbasis, and the operator fx, A G [0, tt), represents a 

rotation by an angle A about the z-axis. Since any such rotation leaves the standard 
basis (IJz)} unchanged, the second MU basis can be transformed into The 
matrix representation of the resulting pair of MU bases reads 

All other pairs of MU bases of the space are, in fact, equivalent to this one. 
A third basis MU to these two bases consists of the states given in Eq. f|TT]) if 
A = ±7r/2, producing Thus, all pairs of MU bases in are equivalent to 

{\jz)', \jx)}, and all triples are equivalent to {\jz)', \jx)', \jy)}, as is well known. 

In dimension three, one of two given MU bases can always be mapped to the standard 
basis {|J^),J = 0,1,2}, so that the second basis consists of states of the form 

\A) = ^{\0,) + e'^lh) + , e,^e[0,27r), (13) 

exploiting the fact that the overall phase of a quantum state has no physical meaning. 
One can construct three states of this form which are pairwise orthogonal: writing 

|A^) = i=(|0,)+7e*«IU+5e^^''|2.)), \^\ = \6\ = 1 , (14) 

the condition (AIA-*-) = implies 7 + 5 = —1. A geometric argument in the complex 
plane implies either 7 = 0; and 5 = w^, or 7 = and 5 = co, where u = e^'^*/^ is a 
third root of unity. We denote the resulting basis by 

{|A),|A^),|A^)} = {4,|J,.)}, (15) 

where the triple {| J^)} = {\Jx), = 0, 1, 2} consists of the eigenstates of the shift 
operator X3, and the operator R^^^ is diagonal in the z-basis such that \A) = i?^ ,,10^.), 
cf. Eq. f|T3|) . The free parameters ,^,77 in the pairs of MU bases {\Jz)] f^£,,rf\Jx)} can 
be removed by a suitable redefinition of the phases of the states in the standard 
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basis Thus, all pairs of MU bases of are equivalent to the pair {\Jz)', \Jx)} 

which may be represented by 



{I;F,}={ \ 1 ];-={ 1 CO cu' ]} , (16) 





where = is the Fourier matrix in C^. Note that two more orthonormal bases 
of states MU to the pair {\Jz)] \Jx)} emerge if one sets either e'^ = e"^"^ = u or 
gj? ^ g«;y = ^2 f[T^ . We will denote these bases by {\Jy)} and {\Jw)}, 

respectively, and their matrix representations are given by 

{)• "-^7s{^: ' '''' 

/ \ 1 / 

which are also MU with respect to each other. The matrices H^, Hy and are 
complex (3 x 3) Hadamard matrices, i.e. they are unitary and the moduli of all their 
entries are equal to 1/ -\/3- 

Two triples of MU bases now result from adding either {| Jy)} or {| to the pair 
{\Jz)'-, \Jx)}- These triples are equivalent to each other as follows from taking the 
complex conjugate (defined in the ^-basis) of the triple {\Jz)'i \ Jx)'i \ Jy)}'- the complex 
conjugation only affects the ordering of states within {| Jx)} while {| Jy)} turns into 
{|Ju,)}. Thus we conclude that the triples are indeed equivalent which we express 
formally by writing 

Consequently, all MU triples are equivalent to the triple {\Jz)'-, \Jx)'-, \Jy)}i and the 
complete set of four MU bases in C'^ is also unique, as is well known. 



3 Constructing product bases in dimensions four 
and six 

The first step towards an exhaustive list of pairs and triples of MU product bases in 
dimension six is to construct all locally inequivalent product bases in C^C^C^. Once 
these are known, the requirement of any two such bases to be MU will impose further 
constraints. It will be helpful to initially carry out this construction in dimension 
four. Thus, we will first derive all inequivalent product bases of the space ® C^, 
followed by a similar construction for a six- dimensional space. 
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3.1 All product bases in d = 4 



We now show that each product basis in d = 4 is equivalent either to the standard 
direct product basis or to a member of two famihes of indirect product bases, each 
depending on two real parameters. Any (orthonormal) product basis in the space 
® must have the form 

\i'2,(p2), \i^3,h), |^4,</)4)}, (19) 

where \ipn), \4'n) £ for n = 1 ... 4. The conditions 

{ipn, 4>n\lpn' , (Pw) = {ipnli^n') {(Pn\(pn') = Snn' , n, = 1 . . . 4 , (20) 

imply that at least two states of the first factor must be orthogonal. However, no 
three orthogonal states exist in C^, so that upon calling {ipi) = \a) we must have 

||a,0i), |a-^,02), |V^3,03), |V^4,04)}, (21) 

with \1lj2) = la-*") being the unique state orthogonal to \a). Now we need to consider 
two separate cases: we can have either \ip-s) = \o) (or, equivalently, \ip3) = la"*")) or 
= \b) such that < < 1, meaning that the state \b) is neither a multiple 

of the state \a) nor orthogonal to it; we call such a vector \b) skew to \a). 

By a simple argument using the restrictions imposed by the orthogonality conditions, 
one finds that three different bases result: 

Bo = {\a,A), \a,A^), \a^,A), |a^,A^)}, (22) 

B^ = [\a,A), \a,A^), \a^,B), |a^,5^)}, (23) 

B2 = {\a,A), \a^,A), \b,A^), (24) 

The basis Bq is a direct product basis while the bases Bi and B2 are not. After 
performing suitable LETs, we can thus summarise the complete list of product bases 
in dimension four as follows. 

Lemma 1. Any orthonormal product basis of the space ® is equivalent to a 
member of one of the families 

2o = {\jz,k;,)}, 

X2 = {|j.,0,),|i)j„l,)}, (25) 

where the operators u,v & SU{2) act on the space such that the states \0z) and 
mIO^), as well as the states \0z) and v\Oz) , are skew. 
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Note that the parameters on which the operators depend have been chosen in such 
a way that no product basis occurs more than once. A number of LETs (cf Sec. 
12. ip have been used to bring the bases into the form given in the lemma. The basis 
Bq in (l22l) has been mapped to Xq by means of a transformation Ui (g) U2 such that 
Ml maps the pair of states {|a), la."*")} to the standard basis {IJz)} of C^, and U2 is 
defined analogously. Thus, the bases Bq and Xq are equivalent to each other. We 
apply a similar transformation to the basis Bi in (123|) mapping two of the bases to 
the standard basis. The freedom to choose a third basis, associated with the pair 
{\B), l-B"*")}, is represented in Xi by the undetermined unitary operator u acting on 
the standard basis. The same reasoning brings B2 into the form fl2S]) except that the 
roles of the two spaces are swapped. Since a complex conjugation reflects points on 
the Bloch sphere about the xz-plane, only half of all the unitaries u (and -0) need to 
be considered in Lemma [TJ In other words, the bases associated with the unitaries u 
and u*, given by the complex conjugate of the matrix representing u in the 2;-basis, 
coincide. 

The symmetry of the space (S) is reflected in the fact that we found two bases 
Xi and X2 which are identical except for the order of the factors. If we stick with 
the idea that LETs dictate whether two product bases are equivalent to each other, 
we need to consider these bases as inequivalent. Thus, the complete set of product 
bases consists of two families each of which depends on two parameters due to the 
5't/(2)-transformations u and v. Not all three parameters of a transformation in 
SU{2) are relevant since the overall phase of quantum states is physically irrelevant: 
each pair of opposite points on the Bloch sphere defines an orthonormal basis of 
so that the set of all bases depends on only two real parameters. Note that the sets 
Xi and X2 of Lemma [1] are both connected to the product basis Xq. 

The symmetry becomes particularly obvious if we represent the bases of Lemma [T]by 
quantum circuits. The idea is to visualise the operation needed to map the states of 
the standard product basis Xq into the desired product basis by means of a quantum 
gate. This is always possible since any two orthonormal bases are connected by a 
unitary operation. Obviously, the trivial gate, described by the identity /, maps 
the four vectors of the standard product basis to itself. Fig. ([1]) shows that (non- 
local) controlled-M and controlled- i) gates are required to output the bases Xi and X2, 
respectively. As expected, the two circuits are identical upon swapping the qubits. 
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u 



Figure 1: Two quantum circuits to create the product bases Ii and I2, respectively; the 
unitaries u and v only act on the target qubit if the control qubit is in the state jl^). 

3.2 All product bases in d = 6 

To construct all product bases in dimension six we use the same method as in 
dimension four. Any product basis in the space ® C'^ takes the form 

|Z^2,^2), |t^3,^3), 1^4,^4), 1^5,^5), 1^6,^6)}, (26) 

with states \tpn) G and I'^n) G for n = 1 . . . 6 , satisfying the orthogonality 
conditions 

{'lpn,^n\'^n','i^n') = {lpn\^n'){'^n\'^n') = 6nn' , Tl, Tl' = 1 . . . 6 . (27) 

The states n = 1...6, in (126|) must contain at least two (not necessarily 

different) pairs of orthogonal states. If they do not, the orthogonality conditions 
require four orthogonal states in C^, which do not exist. In fact, the remaining two 
states in must also be orthogonal, which implies that the product bases of 
will come in three flavours. The states {ipn) ,n = 1 . . .6, fall into three pairs of states 
consisting of either three, two, or only one pair of orthonormal bases. The following 
lemma summarises the results of the detailed arguments given in Appendix [Al 

Lemma 2. Any orthonormal product basis of the space ® C'^ is equivalent to a 
member of one of the families 

Xo = {Ij^, Jz)} , 

Xi = {|0„J,),|l„f/J,)}, 

X2 = 0,), |mO„ 1,), |mO„ 2,), \uU, VU), \ul,, V2,)} , 

T3 = {\jz,0z),\vj.,h),\^Jz,2,)}, (28) 

with j = 0, 1 and J = 0, 1, 2; the operators u,v,w E SU{2) and U,V E SU{3) act on 
and C^, respectively, with V leaving the the state \0z) invariant; the parameters 
of the operators it, . . . ,V are chosen in such a way that no product basis occurs more 
than once. 

Without any restrictions on the five unitary operators u, . . . ,V some product bases 
would occur more than once in this list. For example, if f/ = J, the basis Xi turns 
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into Xo; similarly, the bases associated with U and U* are identical. We could remove 
such multiple occurrences by appropriately restricting the unitary operators but it 
is rather cumbersome to do so and not particularly informative. 

Compared to dimension four, the number of families of indirect product bases have 
increased, and they contain transformations generated by elements of the group 
SU{3). Clearly, there is no scope for symmetry under exchanging the two spaces 
of the product (g) C^. The families Xi to X3 each depend on a number of free 
parameters: Xi has six free parameters due to the unitary U ; two free parameters are 
associated with each SU (2)-transformation present in X3, while X2 is a five-parameter 
family - the transformations due to V, which is effectively an S'f/(2)-transformation, 
brings not only two but three parameters because the overall phase of the states in 
the two-dimensional subspace spanned by and 12^) does not drop out. Figs. 
([2]) and ([3]) show quantum circuits to generate the inequivalent product bases in 
dimension six. 



— u — 



Xi 



Xn 



u 



— V 



X 



xt - 



Xo 



Figure 2: Quantum circuits for a qubit (upper wire) and qutrit (lower wire) to create 
the bases Ii and I2, respectively; the controlled-?/ and controlled-!^ gate act on the qutrit 
only if the control qubit is in the state \lz)'-, the unitary u\ the adjoint of n, acts on the 
qubit only when the control qutrit is in the state \lz)', and the operator X acts as a shift 
on the standard basis of C'^. 
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X 



Figure 3: A quantum circuit for a qubit (upper wire) and qutrit (lower wire) to create 
the basis X3; the unitaries v and w act on the qubit only if the control qutrit is in the state 
and the operator X acts as a shift on the standard basis of C^. 



4 Adding MU product states to sets of orthogonal 
product vectors 

In this section we derive a theorem which will play a crucial role in the construction 
of all pairs and triples of MU product bases in dimension four and six. This theorem 
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is inspired by a constraint on two direct product bases to be MU, obtained in |10| : 

Lemma. Two [direct] product bases {\ja,Ja)} o-nd {|A;b, i^b)} in dimension d = pq 
are MU if and only if \ja) is MU to {kh) in dimension p and | J^) is MU to \Kb) in 
dimension q. 

This result does not cover indirect bases. To (partly) remedy this shortcoming, we 
will present two different ways to generalise this Lemma. Firstly, we find a constraint 
on each product vector if it is to be MU to a specific set of product vectors; this 
result is obtained for spaces of arbitrary composite dimension d = pq. Secondly, we 
derive constraints on a product vector required to be MU to any (direct or indirect) 
given product basis of the spaces C'^ or C^. 

Consider p product states i = 1 . . . with an orthonormal basis {l-i/^j),-? = 

1 . . .p} of the space C^, and with G C''. After swapping the two factors in Eq. 

the product basis {\jz,Oz), \ jx, ^z)}, ioi example, is seen to consist of two sets 
of this form. We find that only particular product states can be MU to such sets of 
product states. 

Lemma 3. The product state |0, $) in dimension d = pq is MU to the set of 
orthogonal product states {{i/ji, '^),i = 1 . . .p} if and only if \(f)) is MU to {tpi) G 
and 1$) is MU to |^) G C. 

If = 1/p and = 1/q, then the product states are indeed MU in the 

space CP'' since it follows that K?/;^, $) p = | p = 1/pq. To prove 
the converse, we assume the product states are MU, | {ipi, $) p = 1 /pq. Summing 
over i = l...p, we obtain |(\E'|$)p = 1/q upon using the completeness relation 
Yli \ = 1- This result immediately implies that = l/p,i = 1 . . .p, 

also holds. 

Note that one can swap the roles of the factors in the tensor product. Then Lemma 
[3] restricts the form of any product state which is MU to a set of q orthogonal states 
= 1 • • • I?} with an orthonormal basis {|\E^j), i = 1 . . . g} of the space O, 
and with {ifj) G C^. 

This result covers the Lemma given at the beginning of this section. To see this, 
group the basis {\ja, Jo)} into q sets of p orthonormal vectors {\ja, la)}, {\ja, 2a)} • • • 
{\jai Qa)}', then, by Lemma [21 any product state |0, $) is mutually unbiased to each 
set of vectors if and only if the state |0) is MU to all states \ja), and the state |$) 
is MU to all states \Ja)- By replacing the state |0, $) with a vector from the basis 
{\ka, Ka)} and repeating the argument for all states in this basis, one arrives at the 
Lemma for direct product bases. 
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The following generalisation uses the fact that we know all direct and indirect prod- 
uct bases in dimensions four and six. 

Theorem 1. The product state |0, $) G C^^d = pq <6, is MU to the product basis 
^fj)} with i = 1. . .pq, if and only if |0) is MU to G and |$) is MU to 

We prove this statement by considering the cases d = 4 and d = Q separately: 

• d = A: All product bases in dimension four are given by the bases Xo,Xi and X2, 
collected in Lemma [1] Each of these bases can be divided into groups of states of the 
form {{ipj, = 1, 2}, or {\ip, ^ = 1, 2}. Thus, Theorem [1] follows immediately 
from Lemma [31 

• d = 6: It is sufficient to consider the four families of bases given in Lemma [2J 
Each of the bases Xq, Xi and X3 can be split into sets of the form required to apply 
Lemma [3l thus, Theorem [1] holds for these bases. To complete the proof, we need 
to consider the basis I2 which has no such decomposition. To begin, suppose that 
the basis X2 is MU to the state |0, $). According to Lemma [3] this state is MU to 
the pair {|j„0,)} if both |(0|O,)p = = 1/2 and |(<f|0,)p = 1/3 hold. The 
state 10, $) also needs to satisfy 

|(0|nO.)n($|UP = |(0|nO.)n($|2.)p = i; (29) 
adding these two constraints we find 

muo,)\'{mh)\' + = 1- (30) 

Using X]j I ('^l^'^z) P ~ ^-^^ completeness relation of the basis {|</^)}, and 

|($|0^)|2 = 1/3, we find that |($|1^)P + |($|2^)p = 2/3. Substituting this identity 
into ([30]) leaves us with \{(f)\uO,)\'^ = 1/2, so that |($|1^)P = \{^\2,)\'^ = 1/3 as 
well. A similar argument applied to the pair {\ulz,Vlz),\ulz,V2z)} shows that 
indeed |((/'|m1^)P = 1/2 and |($|1/1,)|2 = |($|\/2^)|2 = I/3, which confirms that 
the state |0, $) is of the desired form. The converse direction of the statement is 
straightforward. 

We conjecture Theorem [1] to hold for all product dimensions d = pq, i.e. d = 
4, 6, 9, 10, . . . However, a proof similar to the one for d = 4,6, would rely on the 
structure of all product bases in composite dimensions d > 6 - which is not known 
to us. 
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5 MU product bases in dimension four 



5.1 All pairs of MU product bases 

To construct pairs of MU product bases in the space C^, we check all possibilities to 
form MU pairs of the product bases displayed in Lemma [H of Sec. 13.11 We find two 
families of locally inequivalent MU product bases given in Proposition [T] below. The 
derivation of Proposition [1] relies on a technique also used for the six- dimensional 
case presented in Appendix [Bl 

Proposition 1. Any pair of MU product bases in the space ® is equivalent 
to a member of the families 

"^0 ^ = kz)] \jx, kx)} , 

^ = Wz, kz), \lz, Sf,kz); \jx, Ox), \rujx, Ix)} , (31) 

where j,k = 0,1, and the unitary operator rotates the basis {\jx)} = {1^^)} = 
{|±)} into the xy-plane according to ri/|±) = (|0^) e^^yiz)) / \/2 for v G (0,7r); the 
operator generates rotations about thex-axis, i.e. s^\kz) = {\0x) + {—l)'''e'^^\lx))/V2 
for n e [0,7r). 

The pair Vq'^^ is the Heisenberg-Weyl pair consisting of two direct product bases. 
The pair of MU bases ^ is a two-parameter family and may contain direct and 
indirect product bases. Notice that the operator can act as the identity since the 
first basis of v[^'' may be the standard basis kz)}- 

The pair v[^'^ turns out to be equivalent under non-local transformations to the 
Fourier basis as follows from mapping the first basis to the standard basis {\jz, kz)}- 
Thus, we have obtained all known pairs of MU bases in dimension four (cf. Sec. 3 
of |14] ) in spite of limiting ourselves initially to MU product bases only. 

5.2 All triples of MU product bases 

Now we are in a position to derive all triples of MU product bases in dimension 
d = 4: we need to determine which of the pairs of MU product bases given in 
Proposition [1] can be extended by a third MU product basis. 

It is easy to see that the MU pair Pq^'' = {\jz,kz)] \jx,kx)} can be extended by 
adjoining a third direct product basis, namely \jy,ky), resulting in the standard 
Heisenberg-Weyl triple. This is the only possibility, as follows immediately from 
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Theorem [D a product state |0, $) is MU to both {\jz, kz)} and {\jx, kx)} only if |0) 
is MU both to {\jz)} and {\jx)}, and if 1$) is MU both to and {\kx)}. 

The pair v['^^ of MU bases cannot be extended, not even by a single MU product 
state. To extend the pair by an MU product state, one would need to find a state in 
which is MU to the three bases {j/i^x)} and {\r^k^)}. Since u G (0,7r), no 

two of these three bases coincide and there is no state in the space simultaneously 
MU to three distinct bases. As a consequence, the number of MU product triples is 
rather limited in dimension four. 

Proposition 2. Any triple of MU product bases in the space (g) is equivalent 
to 

"Tq ^ — {1^2! ^2)) \jx^kx)', \jyjky)} . (32) 

Using Theorem [1] again, the non-existence of even a single product state MU to the 
triple Tq'^^ follows immediately — all states MU to the triple must be entangled. 

This observation agrees with results reported earlier. For the two-qubit system 
considered here, a construction of the five MU bases based on the Galois field GF{4) 
has been given in [18]. The complete sets obtained turn out to be equivalent under 
local unitary transformations, and they necessarily consist of three MU bases made 
up from separable (i.e. product) bases while the remaining two contain maximally 
entangled states only. This structure also emerges from an approach which exploits 
the fact that any complete set of MU bases of a bipartite system in C'' contains a 
fixed d-dependent amount of entanglement [10]. When d = A, this result implies 
that for a complete set of MU bases containing the triple Tq^^ , the other two bases 
of the quintuple must consist of entangled states - in fact, only maximally entangled 
states are permitted. In [11], the entanglement structure of complete sets of MU 
bases related to Heisenberg-Weyl operators in prime-power dimensions has been 
studied leading to a generalization of the result for dimension d = 4: in bipartite 
systems of dimension d = p^ a. number of (p + 1) MU bases must consist of product 
states while the remaining ones contain only maximally entangled states. 

6 MU product bases in dimension six 
6.1 All pairs of MU product bases 

We will now construct all pairs of MU product bases in dimension six following the 
method used in dimension four (cf. Sec. 13. ip . To obtain a MU pair we take each 
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basis listed in Lemma [2] and go through all possibilities of adding one of the product 
bases Bo to B3 (cf. Eqs. fl471l46]l50]l5T]) of Appendix . 

When constructing pairs of MU product bases, it is not necessary to include the 
basis X2 in Lemma [2l We will show now that the operator V must either act as the 
identity on the pair of states {jlz), 12^)} or swap them, i.e. only a = or /3 = 
are allowed in the expression V\lz) = <y\lz) + f^l'^z)- However, in both cases the 
simplified product basis X2 turns into a special case of X3, given in (128|) . 

Here is the reason why the operator V must simplify in the way just described. 
Apply Theorem [1] to the product state !</),$) required to be MU to I2: the state 
1$) must be MU to all six vectors of present in X2. Consequently, all states in 
C'^ which occur in the bases Bq to B^, defined in Eqs. fl47ll46|50|l5T]) - these are all 
candidates for a second product basis MU to X2 - must be MU to the standard basis 
{| J^)} of C'^. Now, each of these four bases contains another orthonormal basis of 
C^, namely {\A), \A-^), \A-^)}. There is a two-parameter family of such states, given 
in Eq. (IT^ . However, these states must also be MU to the state V\lz) of the basis 
X2. For the states \A) and \A^), this requirement reads 

\{A\{a\lz) + = \{A^\{a\lz) + m)\' = ^- (33) 

Now using the explicit expressions of the states \A) and 1^-*-) given in Eq. (fT6|) and 
the identity |ap + |/3p = l, the first equality leads to 

\l + u\\a\\(3\ = \a\\(3\, (34) 

which implies that either a = or /3 = 0. Thus, for the construction of pairs it is 
sufficient to use the restricted basis 

A = {\jz,Oz), \ujzAz), \u3z,2,)} (35) 

instead of X2 given in Lemma [21 All bases of this form, however, are contained in 
X3 if one chooses v = w = u in fl28|) . This simplification also holds for the basis B2 
when occurring in a pair of product bases. 

The actual derivation of all MU product bases in dimension six is lengthy but 
straightforward. The calculations have been relegated to Appendix [B] except for 
the pairing of the basis Xi with Bi, which gives rise to the pair V3. The proof that 
no other (non-trivial) pair of MU product bases results from {Xi; Bi} has been ob- 
tained by A. Sudbery, and it is given in Appendix O We now summarise the results 
derived in these two appendices. 
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Theorem 2. Any pair of MU product bases in the space ® C'^ is equivalent to a 
member of the families 

= {\jz, Jz); \jx, Jx)} , 

T^l '^z}i |0x; Jx) ^ -^^,r)Jx)"\ ) 

"P2 = {|0z, ^z), II2, Jy)] |0^., Jx), jlx; Jw)} ! 

■p3 = {10^, Jz), |lz, 'S't^^^^Jz); Ox-), |r^jx, Ix), 2^)} , (36) 

with j = 0,1 and J = 0,1,2. The unitary operator R^rj is defined as R^^i = 
|02)(02| + e*^|l2)(lz| + e*''|22)(22| , /or ?7,^ G [0,27r), and Sc^^^ is defined analogously 
with respect to the x-basis; the unitary operators fcr and fr act on the basis {\jx)} = 
{|±)} according to f,|j,.) = {{0^) ± e'''\U))/V2 for a E (0,7r), etc. 

As before, the ranges of the parameters are assumed to be such that no MU product 
pair occurs more than once in the hst. The pairs Vq and V2 have no parameter 
dependence, the pair Vi depends on two parameters, while V3 is a four-parameter 
family. 

Theorem [2] represents the first main result of this paper. It states that there are 
continuously many possibilities to select pairs of MU bases which, however, can be 
listed exhaustively. In the remainder of this paper we will proceed by analytically 
constructing all triples of MU bases which exist in = 6. This will lead to our 
most important result, namely Theorem 4 in Sec. [7] which states the impossibility 
to extend any MU product triple by even a single MU vector. Thus, complete sets 
of MU bases in d = 6 will contain at most pairs of MU product bases. 

An alternative method to exploit Theorem [2] has been pursued in [20]. Upon using 
suitable non-local unitary transformations and known results obtained by computer- 
algebraic methods, the strongest possible statement about MU product bases is then 
derived: if a complete set of seven MU bases exists, it will contain at most one 
product basis - which may be chosen to be the standard basis. 

6.2 All triples of MU product bases 

It is straightforward to enlarge the existing pairs of MU product bases in Theorem [2] 
to triples: simply add the MU product bases listed in Lemma [2], one after the other, 
to each of the pairs Vq to V3 and check whether a valid MU product triple results. 

Neither of the pairs V2 and V3 in Theorem [2] can be extended by a single MU product 
state. To do so, we would need a vector MU to the three distinct bases {|j2)}j 
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and {{"Tajx)} in the space C^, or a vector mutually unbiased to four MU bases in the 
space C^. No such states exist, implying that any state mutually unbiased to these 
pairs must be entangled. 

The pairs Vo and Vi can be extended by a further MU product basis since there 
exist vectors of the spaces and that satisfy the necessary conditions. To 
obtain the complete list of all MU product triples in we thus need to search for 
possible extensions of these two pairs by a third product basis. Starting with Vq, it 
is possible to extend this pair by either Bo or Bi. 

• {Vo]Bo}: If we choose the third basis to be of the form Bq, there are only two 
choices, {\jy,Jy)} or {\jy,Jw)}- Using the local complex conjugation I ® K, the 
resulting triples are found to be equivalent, 

{\3z, Jz)\ bx, Jx)\ by, Jy)} ~ {|jz, Jz)] \jx, Jx)\ Jw)] ] (37) 

consequently, all triples of this type are equivalent to the Heisenberg-Weyl triple 

To = Jz)] IJx, Jx)] by, Jy)] ■ (38) 

• {Po; Bi}: If we extend Vq by an indirect product basis of the form Si, there are 
only two choices, {|0y, Jy), J^,)} or {|0y, J„,), J^)}. Again, a local complex 
conjugation I maps one of the triples into the other, 

{b'^, Jz)\ \jx, Jx); \0y, Jy), \ly, J^)} ~ { b ^ , Jz)] \jx, J x) ] \0y, ^), Jy)} , (39) 

leaving us with the triple 

Tl = Jz); bx, Jx); |0y, Jy), \ly, J^u)} ■ (40) 

Now turning to the pair Vi, we again attempt to obtain a triple by adding either 
Bo or Bi. 

• {Vi, Bq} or {Vi, Bi}: First, extend the pair Vi by a direct product basis, resulting 
in either {b^, Jz); \0x, Jx), Jy); \jy, Jw)} or J^); \0x, Jx), \lx, Jw); \jy, Jy)}- It 
is not difficult to apply suitable LETs to transform them into the triple 7i. Now 
extend the pair Vi by an indirect product basis Bi. This leads to a contradiction 
since we would need the states {|i?^^Jj,)} in Vi to coincide with which is 
not allowed. 

This completes the construction of all MU product triples in dimension six, leading 
to the second main result of this paper. 
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Theorem 3. Any triple of MU product bases in the space (g> is equivalent to 
either 

To = {\jz^ Jz); |jx, Jx); \jy, Jy)} , 

or Tl = {\jz, Jz); Jx); |0y, Jy), \ly, J^)} • (41) 

According to Theorem [H neither of these triples can be extended by a single MU 
product state. Thus, any complete set of seven MU bases in dimension six will 
contain at most three product bases, and if it does, the triple must be equivalent to 
one of those in Theorem [3l In the following section we will obtain an even stronger 
result. 

7 Excluding triples of MU product bases from com- 
plete sets 

In this section we derive the third main result of this paper. 

Theorem 4. No triple of MU product bases in dimension six can be extended by a 
single MU vector. 

In other words, no complete set of seven MU bases in = 6 contains a triple of MU 
product bases. This result relies on a computer-algebraic proof in jS] , which finds a 
total of 48 vectors MU to the pair of eigenbases of the Heisenberg-Weyl operators 
Xq and Zq, giving rise to sixteen different orthonormal bases. However, none of 
these bases allows one to extend the given pair beyond a triple of MU bases. 

The present construction of MU product triples effectively produces twelve (and only 
twelve) product vectors that are MU to the pair Vq = {\jz,Jz); \jx,Jx)}, namely 
{\jyi Jy)} ctiid {\jy, Jw)}, and they give rise to the only two inequivalent triples of 
MU bases. To and 7i. Since Vq is equivalent to the eigenbases of Xq and Zq, clearly 
these twelve product vectors must figure among the 48 vectors given in [6]. 

To show this, we must first deal with a difference in our definition of the HW 
operators. The HW pair used in jB] does not have the same form as Vq since the x- 
basis in [B] is the eigenbasis of the operator Xq, whereas we have used the eigenbasis 
of the operator X2 ® X3 (cf. Eq. ([5])). Nevertheless, both pairs of bases turn out 
to be equivalent using a non-local unitary transformation. By writing the operators 
as matrices, we find that X2 <S) X^ = P25XQP25, where P25 is a permutation matrix 
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permuting rows two and five. This non-local transformation brings the eigenbasis 
of Xq into product form, i.e. {\jx, Jx)}, by multiplying it with P25 from the left. 

The same transformation must also be applied to the list of 48 vectors so that 
they are MU to the pair Vo- After multiplying each of these vectors by the matrix 
P25 from the left, one easily identifies the twelve product vectors, numbered by 
1, 2, 5, 6, 9, 10, 13, 14, 17, 18, 21 and 22 in the Appendix of the updated version of [6]. 
For example, the vector labelled (1) transforms as follows: 



—P25{1, a^, 1, -a^, -o?, -o?Y =—={\, -o?, 1, -o?, a^, -o?Y 
v6 v6 



^ — {l^-^Y®{\^u\\Y (42) 

where a = e^'^*/-^^ and to = e^'^*/^. This vector is the product state \ly, ly). 

The twelve vectors give rise to four of the sixteen orthonormal bases which are 
MU to the original pair. These product bases are covered by the product bases we 
construct when extending the Heisenberg-Weyl pair Vq to a triple; however, only 
two of the four triples are locally inequivalent as follows from exploiting suitable 
local equivalence transformations. 

Upon combining the computer-aided result just described with Theorem [3] it is 
straightforward to arrive at a result that excludes product triples from being part 
of a complete set of seven MU bases. The triples of MU product bases % and 7i 
both contain the Heisenberg-Weyl pair Vq = {\jz, Jz)] lix, Jx)}^ and it is impossible 
to extend this pair by more than a single MU basis according to |6]. Since Theorem 
[3] provides an exhaustive list of MU triples in the space C^, it follows that no 
complete set of seven MU bases in d = 6 contains a triple of MU product bases. 



8 Summary and discussion 

By limiting ourselves to orthonormal product bases, we have been able to obtain a 
number of analytic results regarding the existence of MU bases of the space ® 
C^. After identifying all orthonormal product bases of this space, presented in 
Lemma [21 we have constructed an exhaustive list of pairs of MU product bases. 
They come in four different flavours according to Theorem [2l Next, Theorem [3] 
states that, in addition to the Heisenberg-Weyl triple, there is only one other locally 



21 



inequivalent triple of MU product bases. The absence of quadruples of MU product 
bases agrees with Zauner's conjecture [19] that there are no more than three MU 
bases in dimension six. 

The derivation of the list of MU product pairs and triples has been simplified con- 
siderably by the content of Theorem [TJ It spells out severe restrictions on the form 
of product states required to be MU to certain sets of orthonormal states in the 
space ® C"^. We have established Theorem [T] for dimensions d = 4 and d = 6 
only, since the proof relies on enumerating all orthonormal product bases in these 
dimensions. 

Theorem [3] allows us to partly replicate results obtained by means of a computer- 
algebraic method. Out of the 48 vectors mutually unbiased to the Heisenberg-Weyl 
pair Vo, found in jS], we successfully recover twelve, and they are shown to be 
equivalent to product vectors. 

The most important consequence of exhaustively enumerating MU product bases in 
dimension six is a bound on their allowed number in complete sets of MU bases. 
Applying Theorem [T] to the triples of MU product bases in (8> C^, namely To 
and 7i, directly implies that no single product state can be MU to any of them. 
However, a stronger result is within reach, spelled out in Theorem |H it is impossible 
to complement either To or 7i by any MU vector. This follows from combining 
Theorem [3] with the results derived in j^. Thus, a complete set of MU bases in 
dimension six cannot contain a product triple. This is in marked contrast to the 
prime-power dimension where a complete set of MU bases necessarily contains 
{p + 1) MU product bases constructed from the tensor products of Heisenberg-Weyl 
operators [11]. The exhaustive list of MU product pairs given in Theorem 2 has 
been used to derive a result even stronger than Theorem 4, reducing the number of 
allowed MU product bases to just a single one |2U| . 

A similar situation has been described in |21| where a different class of MU bases 
is studied. Given a "nice unitary error basis", consisting of d"^ suitable matrices, 
one can search for MU bases within these sets. In the case of dimension six, it is 
shown that any partition of a nice error basis gives rise to no more than three MU 
bases. This limitation and the non-existence of more that three MU product bases 
are independent results: MU product bases and MU bases arising from nice unitary 
error bases are structurally different. For example, our construction reproduces the 
continuous family pj^'' of MU product pairs in ci = 4, and it is known that some of 
the pairs in this family are inequivalent to MU bases stemming from nice unitary 
error bases |22| . 

Our considerations are backed by deriving corresponding results in the Hilbert space 
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of two qubits, i.e. C^. In this case there exists a symmetry between the 
two factors and the enumeration of MU product pairs and triples is much simpler. 
Clearly, when a qubit is combined with a qutrit, no such symmetry exists. We 
believe that the symmetries between the subsystems present in only prime-power 
dimensions are the ultimate reason that additional "identities" exist which allow for 
the construction of complete sets of MU bases. 

Let us conclude by formulating a conjecture which emerges naturally from our re- 
sults: we expect Theorem [1] to hold for all composite dimensions d = pq > 4, not 
only for c? = 4 and d = 6. Our pedestrian proof in these dimensions relies on 
enumerating all orthonormal product bases. However, the set of product bases in 
composite dimensions is likely to possess a certain structure which, once spelled out, 
should allow for a more elegant proof applicable to arbitrary composite dimensions. 
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A Appendix 

In this Appendix we derive all product bases in dimension six reported in Lemma [2l 
The six states iV'j) = 1 • • • 6, of a product basis {It/^j, defined in (12^ must 

contain at least two (possibly identical) pairs of orthogonal states. If there was 
only one pair (with the remaining four states of the space non-orthogonal), the 
orthogonality conditions f l27|) would require four orthogonal states \^ j) G C'^, which 
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do not exist. Thus, denoting the orthogonal pairs by {|a), \a^)} and {\b), l^"*")}, the 
product basis must take the form 

{\a,^,), \a^,^2), \b,^3), \b^,^4), 1^6,^6)}. (43) 

The states {ipr,) and must also be an orthogonal pair. To see this, assume 
that they are skew (or identical) and they are both skew to the states \a) and 
1 6); then the state l^&e), for example, must be orthogonal to the orthonormal triple 
{|\E'i), I^E^s), 1^5)}, which is impossible. Here we have assumed that \a) and \b) are 
not orthogonal; if they are, we use the orthonormal triple |^4), l^s)} instead. 

The same conclusion can be drawn if the states and l^/'g) are skew (or identical) 
but one of them coincides with any of the four states |a), \b), \a-^) or l^-*-). Thus, we 
conclude that any product basis of the space must be of the form 

{\a,^,), \a^,^2), \b,^3), 16^,^4), |c,^5), |c^,^6)}. (44) 

Now it is obvious that we need to consider three different possibilities depending on 
how many of the three bases of the space coincide. 

Case 1: If all three bases coincide, we have 

||a,^i), |a^,^2), |a,^3), |a^,^4), |a,^5), |a^,^6)}- (45) 

These six states are orthogonal only if the three states l^l/s) and |\l/5) are or- 

thogonal to each other, as well as the triple {|^l/2); 1^4); l^e)}- Upon denoting the 
first triple by {\A), \A-^), \A-^)}, where \A-^) G is a vector orthogonal to \A) and 
lA"*"), we obtain 

Bi = {|a, A), |a, A^), |a, A^), |a^,5), |a^,5^), |a^,5^)} , (46) 

also introducing an arbitrary second triple of orthogonal states. If the two triples 
coincide, we find the important special case of a direct product basis 

i3o = ||a,A), \a,A^), \a,A^), \a^,A), \a^,A^), |a^,A^)}, (47) 

Case 2: If only two of the bases in coincide, we find 

{\a,^,), \a^,^2), \b,^3), |&^,^4), |&,^5), |&^,^6)}. (48) 

As in Case 1, each of the triples l^a), l^s)} and {|\I^2), |^1'4), l^e)} must be an 

orthonormal basis of C^. However, we also need to have 

(^i|^4) = (^il^e) = 0, (49) 
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which means that |\l/3) and |\l/5) span the same subspace as |\l/4) and l^l/e)- It follows 
that l^'i) = |^'2); upon calling this state \A), we are led to a new class of product 
bases of given by 

B2 = {\a,A), |a^,A), \b,A^), \b^,VA^), \b,A^), (50) 

where V^l^-^) = a\A-^) + (3\A-^) and V\A-^) = l3\A-^) - a\A-^), i.e. V is any unitary 
transformation of the two-dimensional subspace of C'^ orthogonal to the state \A). 

Case 3: Finally, we consider the case where the three bases of the space present 
in (jH]) are all different, meaning that |a), \b) , and |c) are pairwise skew. Then, 
the orthogonality conditions directly imply that the triples (I^E'i), l^s), l^s)} and 
{ 1^^2)5 1^4)) I^g)} of orthogonal states must coincide. This leaves us with bases of 
the form 

B, = [\a,A), \a^,A), \b,A^), \b^,A^), |c,A^), |c^,A^)}. (51) 

These three cases complete the construction of all product bases in dimension six. 
Using local equivalence transformations in analogy to the procedure used in Sec. 
13. H one can write the four sets of product bases as displayed in Lemma [21 

B Appendix 

In this Appendix we derive all pairs of MU product bases in dimension six by pairwise 
combining the ortho normal product bases Bq to B3, defined in Eqs. fl47ll46|50|l5T]) . 
In principle, we need to look at only 10 of the 16 pairs {Bf, Bj}, z, j = . . . 3, since 
the order of the bases does not matter: the pairs {Bf, Bj} and {Bj] Bi} are equivalent 
for i ^ j. Using local equivalence transformations, each pair can be brought to the 
form {li] Bj}, i < j, where the bases Iq to X3 are those listed in Lemma [2l As shown 
in the main text, it is not actually necessary to consider the bases X2 and B2 at all, 
reducing the number of cases to six. Parameter ranges are assumed so that no pair 
occurs more than once. 

• {Xq; Bq}: First we extend Xq to a pair of MU bases by combining it with 

Bo = ^^\a,A), \a,A^), \a,A^), \a^,A), \a^,A^), |a^,A^)}. (52) 

The states of Bq are MU to those of the basis Xq if the pair of states {|a), la"*")} 
is any basis of associated with opposite points on the Bloch sphere, i.e. \a) = 
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(|0z) + e^'-'\lz)) / etc., and if the orthonormal basis {\A), is defined as 

in Eqs. ([HD of Sec. EJl 

A local transformation allows us to rotate the states {|a), \a^)} into and to 

simultaneously change the basis {\A), \A^), \A^)} into the basis {\J^)} of so that 
we end up with the known Heisenberg-Weyl MU pair of direct product bases, 

Po = {|j.,J.); (53) 

• {Xo;i3i} and {Xq]B^}: These cases will be covered by the pairs and 
{Xi; i^a}, respectively, since we can treat the basis Xq as a subset of Xi. 

We now construct the three pairs of indirect product bases that contain Xi = 
{\^zi Jz) A'^ziU Jz)} as the first basis, where the unitary tj maps the basis {jJ^)} 
of the space C'^ to another basis. 

• {Xi; Bi}: In a first step, we act with a local unitary on the second basis 

Bi = ||a,A), |a,A^), |a,A^), \a^,B), \a^,B^), |a^,fi^)} (54) 

to rotate the a-basis of states that are MU to {Ij^)} into the basis while the 

74-basis turns into {| Ja;)}, as before. This maps Bi to 

{|0,, J,), |1,,?7'J,)}, (55) 

where we have introduced a unitary U' which parameterises all orthonormal bases 
of relative to the x-basis. The requirement that the states of the pair {Xi; Bi} 
be MU now turns into the problem of identifying all pairs of orthonormal bases of 
C'^, namely {\Jz)] \UJz)} and {\Jx)', \U'Jx)}, such that all states of one set are MU 
to those of the other, viz. 

\{Jz\U'J.)\^ = \{UJz\J.)\^ = \{UJz\U'J.)\^ = I , (56) 

while \{Jz\Jx)\'^ = 1/3 holds by construction. It is easy to see that these conditions 
are satisfied if the bases in (at least) one pair coincide or all four are different, i.e. 
they use up a complete set of MU bases in C'^. In Appendix [C] we present a proof, 
due to A. Sudbery, that these are the only solutions of the constraints f l56l) . Thus, 
if Xi is the standard basis {\jz, Jz)}, then we obtain the MU product pair 

Vi = {\jz,Jz); |0.,J.),|1..,4,J..)}, (57) 

with {\R^,riJx)} defined in Eq. ( ITSl) . However, if we use the complete set of MU 
bases in we obtain the MU product pair 

V2 = {|0„ Jz), \lz, Jy); \0x, Jx), |lx, ^)} . (58) 
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• {Xi; B3}: The second basis reads explicitly 

Bs = [\a,A), \a^,A), \b,A^), \b^,A^), |c,A^), |c^,A^)}, (59) 
and suitable LETs map it to 

{|jx,Ox), |f^jx, Ix), \frjx,2^)} , (60) 

which involve two rotations of the basis about the 2-axis, -r^ and f^-. The 

operator U in Xi must be chosen such that {\UJz)} is MU to the x-basis. All such 
f/(3)-rotations are given by the two-parameter family 

Scx = |0..)(0,| +e^'^|U(l,.| + e*^|2,.)(2,.| , (61) 

diagonal in the x-basis, and defined in analogy to R^ ,^ in Eq. (fT5|) . Altogether, we 
obtain a four-parameter family of MU product pairs, 

^3 = {|0z, Jz), \lz, S(^,xJz)] \jx, 0^), \fjx, 1^.), 2^.)} . (62) 

• {X3;i33}: No pair results when we combine the product basis B3 with X3. The 
standard transformations to simplify B3 lead to 

{|jx,Ox), \rajxAx), |^rjx,2x)}, (63) 

since both the 6-basis and the c-basis must be MU to the standard basis. The only 
basis MU to the three bases {|jx)}; il^ajx)}, and {Irrjx)}, is the standard basis 
which is also true for the case {\rajx)} = {If'rjx)}- Consequently, this would 
force the operators u and v to be the identity, in contradiction to the assumption 
that the three bases of present in X3 do not coincide. 



C Appendix 

Here we report a proof by A. Sudbery that the conditions of Eq. f l56l) in Appendix 
[Bl are only satisfied if the bases in (at least) one pair coincide or all four bases are 
mutually unbiased. If Bi and B2 are orthonormal bases, we write Bi fi B2 to mean 
"Si and B2 are mutually unbiased". 

Theorem 5. Suppose Bq, Bi, B2, B^ are orthonormal bases of satisfying 

{Bo,B,} ^ {B2,B^}. 

Then either Bq and Bi are equivalent bases or B2 and S3 are equivalent bases or all 
four bases are mutually unbiased. 
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Let Bo, Bi, B2, Bs be represented by unitary matrices J, U, V, W, respectively, where 
we have chosen Bq to be the standard basis of C^. We regard the bases U, UP and 
UD, with P a permutation matrix and D a diagonal, as equivalent bases. Note 
that if two ort honor mal bases in C"^, represented by unitary matrices U and V, are 
mutually unbiased, then WV (where the dagger denotes hermitian conjugation) is 
a complex Hadamard matrix H. We can write any (3 x 3) Hadamard matrix as 

H = DFD' or DF^D' (64) 

where D and D' are diagonal and F = F3 is the Fourier matrix defined in Eq. (flGll . 

The condition B2 ^ Bq implies the unitary V is a Hadamard matrix, and since 
FP, the basis B2 is equivalent to a basis represented by = DF. Similarly, B3 is 
equivalent to a basis represented hj W = D'F where D' is diagonal. Now 

B2^iBi =^ V^U = KF'^^^L, (65) 

B^fiBi =^ W^U = K'F^^^L' (66) 
where K, L, K' and V are diagonal and F^*) is either F 01 F^ (z = 1, 2). Hence 

U = DFKF^^^L = D'FK'F^'^^L'. (67) 

We will now examine the relationship between U and the diagonal matrices D,K,L 
in the two cases U = DFKFL and U = DFKF^L, respectively. We can assume 
the leading entries of D and L to be du = In = 1 by absorbing two phase factors 
in the diagonal matrix K. 

Lemma 4. Suppose U = DFKFL where D, K, L are diagonal unitary matrices 
with D = diag{l,a, P). Then either U = PE where P is a permutation matrix and 
E is diagonal, or the matrix elements of U are all non-zero and satisfy 

^12^23^1 = ^13^21^2 = ^11^^22%3, (68) 

and a and (5 are given by 

3 U21U22U23 

a = , (69) 

U11U12U13 

P = a^^. (70) 

^^21^i22 

Let K = diag(7, 6, e) and L = diag(l, (, rj). Then 

/a Vc\ 
U = \ ab a(c a-qa I (71) 
\/3c PCa PvbJ 
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where 



a 



Ul + S + e), 



3 



6 = i(7 + a;(5 + a;2e), (72) 
c= l{j + u^5 + ue). 

Suppose one of a,b,c were zero, say a = 0. Then, since 7, 5, e all have modulus 1, 
they must form an equilateral triangle in the complex plane, so either 6 = U'~f and 
e = u'^'y, when 6 = and c = 7, or 5 = u^'j and e = u'j, when 6 = 7 and c = 0. In 
both cases U is of the form PE. 

If none of a, 6, c are zero, then all the matrix elements of U are non-zero and equations 
(EHD, (EHD and (CD]) follow immediately from (ITB . 

Exactly similar arguments prove 

Lemma 5. Suppose U = DFKF^L where D, K, L are as in Lemma^ Then either 
U = PE where P is a permutation matrix and E is diagonal, or the matrix elements 
of U are all non-zero and satisfy 

M11M23M32 = M12M21M33 = M13M22M31, (73) 
while a is given by (169|) and (5 by 

/3 = a^^. (74) 
M21M23 

We now return to eq. fl67j) and consider the four possibilities for {F^^\ F^"^^). 
Case 1: U = DFKFL = D'FK'FL'. 

Let D = diag(l, a, /3), D' = diag(l, a', /?'). Then, by Lemma HI either U is of the 
form PE (when the bases Bq and Bi are equivalent), or 

a^ = a'^ and ^=(^)\ (75) 

/3 \aj 



Hence a' = a or ua or u'^a, so 



















a 





1 °^ ( 














1 

I or I | . (76) 






30 





This gives 

V = DF= \ a ooa oo^a \ , (77) 

W = D'F = ( a ua u'^a I or j ua u'^a a I or j u'^a a ua | . (78) 

In each case the columns of W are a permutation of those of V . Thus either the 
bases Bq and Bi are equivalent or B2 and B^ are equivalent. 

Case 2: U = DFKFL = D'FK'F^L'. 

Suppose U is not of the form PE. Then both Lemmas H] and [5] apply, and U has 
non-zero matrix elements satisfying (l68|l and ( !73|) . As in case 1, let D = diag(l, a, (3) 
and D' = diag(l, a', /?'). Now a and /3 are given by Lemma HI but a' and (3' are 
given by Lemma O Once again we have = a'^, but now f3'/(3 is not determined 
solely by a' /a: 



(3' f a'Y ursU22 



(3 \(y J M12M23 



(79) 



Using (168|) and (173|). 



\U12U23J \U12J \U23, 

_ Ul3 U23U31 ^21^33 U22 U12U31 UuUs2 
Ul2 U21U32 U22U31 U23 U11U33 U13U31 

= 1. 



?0) 



Hence a' /a and /3'//3 are both cube roots of 1. Write a' = (pa, [3' = xP- If X = '/'^ 
then, as shown in Case 1, the columns of V and W are the same, up to permutation, 
and the bases B2 and B3 are equivalent. \i (f^ then two of 1, x, are equal and 
the third is different. The same is true of the sets {l^ujXi^'^'P} and {1, tu^X; '^</'}- 
Hence the sums a = l + x + 0, 6 = 1 + uox + and c = 1 + u'^x + ^'P all have the 
same modulus. For xi^ ■, the product 



V'^W = F^D^D'F 



is a Hadamard matrix and hence the bases B2 and B3 are mutually unbiased. Thus 
in this case, B2 and B3 are either equivalent or mutually unbiased. 
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Case 3: U = DFKF^L = D'FK'FV. 



This is the same as Case 2 with V and W interchanged. 
Case A: U = DFKF^L = D'FK'F^L'. 

This is similar to Case 1, using Lemma |5] instead of Lemma HI The conclusion is 
the same. 

We have now shown that in every case, either B2 and B3 are equivalent or Bq and 
Bi are equivalent or B2 and B^ are mutually unbiased. But the assumptions of the 
theorem are symmetric between the pairs {Bo,Bi} and {^2,-63}, so we can also 
prove that if B2 is not equivalent to B3 and Bo is not equivalent to Bi, then Bq and 
Bi are mutually unbiased and therefore all four bases are mutually unbiased. 



32 



